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ABSTRACT 

When a white dwarf (WD) is weakly magnetized and its accretion disk is thin, accreted material first reaches 
the WD's surface at its equator. This matter slows its orbit as it comes into co-rotation with the WD, dissipating 
kinetic energy into thermal energy and creating a hot band of freshly accreted material around the equator. 
Radiating in the extreme ultraviolet and soft X-rays, this material moves toward the pole as new material piles 
behind it, eventually becoming part of the WD once it has a temperature and rotational velocity comparable with 
the surface. We present a set of solutions which describe the properties of this "spreading layer" in the steady 
state limit based on the conservation equations derived by Inogamov & Sunyaev (1999) for accreting neutron 
stars. Our analysis and subsequent solutions show that the case of WDs is qualitatively different. We investigate 
example solutions of the spreading layer for a WD of mass M = 0.6M Q and radius R = 9 x 10 8 cm. These 
solutions show that the spreading layer typically extends to an angle of #sl ~ 0.01-0.1 (with respect to the 
equator), depending on accretion rate and the magnitude of the viscosity. At low accretion rates, M < 10 18 g s -1 , 
the amount of spreading is negligible and most of the dissipated energy is radiated back into the accretion disk. 
When the accretion rate is high, such as in dwarf novae, the material may spread to latitudes high enough to 
be visible above the accretion disk. The effective temperature of the spreading layer is ~ (2-5) x 10 5 K with 
approximately T e ff oc M 9 ' 80 . This power-law dependence on M is weaker than for a fixed radiating area and may 
help explain extreme ultraviolet observations during dwarf novae. We speculate about other high accretion rate 
systems (M > 10 18 g s" 1 ) which may show evidence for a spreading layer, including symbiotic binaries and 
supersoft sources. 

Subject headings: accretion, accretion disks — novae, cataclysmic variables — stars: dwarf novae — white 
dwarf — X-rays: stars 



1. INTRODUCTION 

In close binaries containing a white dwarf (WD), the com- 
panion star can fill its Roche lobe either through stellar evo- 
lution or angular momentum loss and begin accreting onto 
the WD to form a cataclysmic variable (CV; Warner 1995). 
Before material can reach the WD it must transport angular 
momentum outward via an accretion disk. Observations and 
theoretical studies show that when the disk can radiate the 
internally dissipated energy (Pringle 1981) it will be thin (i.e. 
the vertical height of the disk is much smaller than the WD ra- 
dius; Shakura & Sunyaev 1973; Meyer & Meyer-Hofmeister 
1982) and rotating at nearly Keplerian velocity. If the WD is 
weakly magnetized, the accreted material reaches the surface 
at the equator and then must pass through a transition region 
to settle into the WD and become part of the star. When the 
WD spin is much less breakup, nearly half of the accretion 
luminosity is released in this region, making it as luminous 
as the accretion disk and likely crucial to understanding the 
luminosity from accreting WDs. 

The transition from the accretion disk to the WD surface 
is typically treated in a boundary layer (BL) model, which 
is simply an extension of the accretion disk with additional 
torques provided by the WD to decelerate disk material to the 
WD rotational velocity. In BL studies the important coordi- 
nate is assumed to be radial and the vertical structure is solved 
by assuming a vertical pressure scale height (e.g. Popham 



& Narayan 1995). The geometry of the BL is of a fattened 
disk near the WD surface. Other BL studies solve the two- 
dimensional problem simultaneously in both radial and lati- 
tudinal directions using numerical techniques (e.g. Kley & 
Hensler 1987; Kley 1989a, 1989b). Kippenhanh & Thomas 
(1978) investigated the properties of "accretion belts" on WDs 
by assuming that the accreted material is marginally stable 
to the Richardson criterion. We now want to understand the 
state of recently accreted material on much shorter timescales 
« 10 4 yr), and assess the effects of a turbulent viscous stress 
on the spreading properties. In Inogamov & Sunyaev (1999; 
hereafter IS99), this problem was approached from a new 
angle to study accreting neutron stars (NSs). Their method 
follows the latitudinal flow of matter on the compact object 
which provides information about the spreading area of hot, 
freshly accreted material which is not captured in BL mod- 
els. We now apply these same methods to the case of WDs 
and call this model the spreading layer (SL), to differentiate 
it from BL studies. We solve the conservation equations de- 
rived by IS99 for fluid confined to the WD surface and find 
that the steady state limit naturally attracts toward solutions 
where local viscous dissipation is nearly balanced by local ra- 
diative cooling. This is qualitatively different than the results 
of IS99 who found that the SLs of NSs show a great deal of 
advection of dissipated energy up to higher latitudes where it 
is radiated. This is mainly due to their high accretion rates 
near the Eddington limit. 
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For a fiducial WD of mass M = 0.6M Q and radius R = 
9 x 10 8 cm, our solutions yield a thin hot band of extent 9sl ~ 
0.01-0.1 with an effective temperature of ~ (2-5) x 10 5 K 
over a range of accretion rates from 10 17 - 10 19 g s" 1 , imply- 
ing that the SL should contribute to accreting WD spectra in 
the extreme ultraviolet (EUV) or soft X-rays. We also find the 
SL's dependence on accretion rate and viscosity. These so- 
lutions show that at low accretion rates (M < 10 18 g s -1 ) the 
SL does not appear above the accretion disk, so that radia- 
tion from the SL may be obscured by the accretion disk or its 
winds. Only when M > 10 18 g s -1 , such as in dwarf novae, 
symbiotic binaries, or supersoft sources, will the SL reach 
high enough latitudes to be directly seen. Perhaps the best 
chance for testing the SL model are dwarf novae in outburst 
and superoutburst. Some observations show a fairly constant 
effective temperature over a large change in luminosity (for 
example SS Cyg; Mauche, Raymond & Mattei 1995) which is 
consistent with a SL. Direct comparisons between our simple 
model and observations can be complicated if the SL emis- 
sion is absorbed and re-radiated by accretion disk winds, as 
is likely for systems which show that the EUV is not eclipsed 
by the secondary whereas the WD is (such as for the edge-on 
OY Car in superoutburst; Mauche & Raymond 2000). Never- 
theless, the SL may be important for calculating the underly- 
ing continuum scattered by these winds. We predict a quick 
timescale for spreading, ?sl ~ 50 s, which means that the typ- 
ical mass of the spreading region (< 10 20 g) is much less than 
the "accretion belt" of Kippenhahn & Thomas (1978). Most 
likely, as angular momentum is tranferred between the SL and 
the underlying WD surface, the underlying surface layers will 
spin up similar to what Kippenhahn & Thomas (1978) calcu- 
late. The short spreading timescale also means that following 
a dwarf nova outburst, when accretion onto the WD has de- 
creased dramatically, the SL will spread over the WD surface 
too quickly to be observed as a separate hot component. 

In §2 we solve for the radial structure of the SL using a 
one-zone model. We derive the differential equations which 
describe the SL in the steady state limit in §3, using IS99 as a 
guide. In §4 we discuss the solutions of the spreading equa- 
tions. We study physically motivated analytic estimates which 
capture the essential features of the numerical integrations. 
These provide power law relations which show how the SL 
changes with accretion rate, viscosity, and the properties of 
the accreting WD. We then present the results of integrating 
the equations numerically to undestand the SL in more detail. 
In §5 we discuss observational tests for the SL and consider 
the current data in relation to our calculations. We conclude 
in §6 with a summary of our findings along with a discussion 
about further work which can be done with this model. 



use an "effective" gravitational constant 



geft = 



GM__v\_v% 
R 1 R R 



(1) 



where vg is the latitudinal velocity, 9 is measured with respect 
to the equator, and we assume that the WD is rotating much 
slower than breakup. Since g e g is independent of radius in the 
thin limit, hydrostatic balance is integrated to give 



P = gesy, 



(2) 



where y is the column depth, defined as dy = -pdz, and z is 
the radial coordinate in this plane parallel layer. 

Radiation pressure may play an important role in the hot SL 
so we define the force per unit mass due to radiation flux at 
the photosphere, 



<?rad = 



KF 



(3) 



where k is the opacity, F is the flux, and c is the speed of 
light. If the one-zone layer becomes optically thin we cannot 
assume the flux provides this force and these approximations 
become invalid. With these definitions g e g = g ra( j defines the 
local Eddington accretion rate. We use an opacity dominated 
by Thomson scattering (k = 0.34 g cm" 2 for a solar compo- 
sition with mean molecular weight per electron of p e « 1.2). 
Free-free absorption may also be an important mechanism for 
radiative transport so we check our assumption of an electron 
scattering dominated opacity once we find the numerical so- 
lutions for the SL in §4. 

We assume that the flux is constant throughout the layer. 
Implicit in this assumption is that the majority of the dissipa- 
tion is at the bottom of the layer where hot, rapidly rotating 
spreading material comes into contact with the cold, slowly 
rotating WD. Radiative diffusion results in 



AacT 5 dT 
3k dy ' 



(4) 



where a is the radiation constant. 

Pressure is related to density and temperature by the equa- 
tion of state, 



pkT aT 4 

P=- — + , 

[im p 3 



(5) 



where lc is Boltzmann's constant and p, is the mean molecular 
weight (/i = 0.62 for a solar composition). Equation (4) is 
integrated and combined with equation (5) to give 



2. RADIAL PROPERTIES OF THE HOT SPREADING LAYER 

Following IS99, we first construct a one-zone model to cap- 
ture the radial extent of the hot, rapidly rotating material in 
the SL, derived in a frame on the surface of a WD. The SL 
is assumed to start near the equator in hydrostatic balance. 
The pressure scale height, h = P/pg ~ 10 7 cm (this will be 
shown in §4), is significantly less than the radius of a typi- 
cal WD so we assume a plane parallel geometry with con- 
stant gravitational acceleration. Initially the azimuthal veloc- 
ity of the spreading material, v^, will be nearly Keplerian, 
~ vk = (GM/R) 1 / 2 , so that gravitational acceleration can 
be significantly decreased by centrifugal effects. We therefore 



3kF 



1/4 



,1/4 



and 



pm p 

P=-j^r(g e ff-gmd)y- 



(6) 



(7) 



We use equations (2), (6), and (7) to describe the layer, where 
>■ is the column depth at the base of the spreading matter. This 
material will initially have a much higher entropy than the 
WD surface when it arrives at the equator. As the material 
moves toward the WD pole it will radiate and slow its rota- 
tional velocity until it becomes part of the WD (Figure 1). 
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Spreading layer with nearly Keplerian 
velocity and higher entropy than 
underlying white dwarf 



Accretion 
flow in disk 



Spreading material merges 
with material below after 
about a minute 



Old accreted material 
co-rotating with white dwarf 
and at low entropy 

FIG. 1 . — A schematic diagram of the spreading layer making contact with the white dwarf. The accretion flow arrives at the white dwarf with nearly Keplerian 
azimuthal velocity. As material piles at the equator, rotational energy is dissipated, heating the layer. The pressure which is created at the equator forces material 
toward the pole, creating a spreading layer. The material continues to travel away from the equator until it loses its azimuthal velocity due to frictional dissipation 
and has an entropy similar to the underlying white dwarf. At this point the spreading material becomes part of the slowly downward advecting material on the 
white dwarf surface. 



3. DERIVATION OF SPREADING EQUATIONS 

Since we are primarily interested in how the properties of 
the accreting material change with latitude, it is easiest to 
consider the equations for conserved fluxes (fluxing in the 
^-direction), integrated radially through the SL. These equa- 
tions are taken in the steady state limit, resulting in a set of 
equations analogous to those derived by IS99. We include the 
majority of the derivation here as it will assist our choice of 
boundary conditions which is different than what IS99 found 
for NSs. In our derivation we primarily use notation consis- 
tent with IS99, but make changes when convenient to do so. 

3.1. Conservation Equations 

We first derive the general, time dependent, conservation 
equations in a frame on the WD. For conservation of mass we 
consider the radially integrated density (i.e. the column depth, 

y), 

d 1 d 

— (2nRycos9) + -— (2nRyv g cos9) = 0, (8) 



dt 

which simplifies to 



R80 



(9) 



Rcos9^ + -^riyvg cos 9) = 0. 
at o9 

Conservation of momentum in the ^-direction is given by 

— (2t: Ryvg cos 9) + - — (2nRyv 2 g cos 9) + 2itRy^- sin 9 
ot R o9 R 

= -2ttRcos9-^-\ [ Pdz] -2ttRtb cos 9. (10) 

R d9 yJo J 

which can be rewritten as, 

d d 2 2 

R cos 9 — (y v e ) + ^ (y v g cos 9) + y v sin 9 



>sw— y J Pdz\-RT0Cos9 1 (11) 



where the last term on the left side is a centrifugal piece which 
comes from the gradient of the spherical coordinate basis vec- 
tors (Landau & Lifshitz 1959), the first term on the right is 
the radially integrated hydrostatic force, and r g is the vis- 
cous stress in the ^-direction. The total viscous stress, r, 
is from the turbulence created as high entropy fluid quickly 
rotates against the lower entropy WD material below. We 
parametrize r in terms of a unitless constant, a, so that 

r = apv 2 , (12) 

where v 2 = v 2 g + v^. In Appendix A we consider lower lim- 
its to a provided by ion or radiative viscosities. These are 
found to be fairly small (a < 10~ 6 ), implying that turbulence 
is probably the dominant cause of friction. To be as general as 
possible we treat a as a free parameter and study its affect on 
the spreading properties for the values a = 10" 4 - 10" 2 . This 
range of values is chosen for the physically realistic solutions 
which result for M - 10 17 - 10 19 g s" 1 . If a < 10" 4 the vis- 
cous stress is too slow to halt the spreading, and the spreading 
flow diverges in numerical integrations. For a > 10~ 2 the dis- 
sipation becomes so large that the scale height puffs up to be 
comparable to the WD radius, making the plane parallel as- 
sumptions invalid. In such a case a more detailed study of the 
radial structure is required than what we provide in §2. We 
discuss the available parameter space for values of a in more 
detail in §4 when we consider the boundary conditions. The 
components of r are 



2 y 6 nr. T 

T 9 = apv —=apv e ^v^ + v z 4> , 



(13) 



and 



T4 > = apv <j> ^v 2 + v 2 4> , (14) 
in the latitudinal and azimuthal directions, respectively. 



The pressure, integrated over a scale height, h, is 



J() P m P \ ac J g e ff-gmdJo 

g e s kTy 



' A dy 



(15) 



5geff-,gradM m / 

where both y and T are assumed to be their respective values 
at the base of the SL, and k is assumed to be constant. We see 
that difficulties arise in the one-zone formulation when the SL 
becomes close to the Eddington limit (g e ff = grad)- This must 
be carefully considered when we solve the resulting differen- 
tial equations in §4. 

Next, there is conservation of momentum in the in- 
direction, 

Q 1^ 

—(2irRyv^ cos 9) + -—(2nRyv 4> vg cos 9) 

-lirRy^- sin 9 = -IttRt^ cos 9. 
R 

Simplifying this expression we find, 



(16) 



d 

Rcos9—(yv q 



d 

) + Qgiyv^ve cos 9) -yv^vg sm9 
= -RtsCOs6, 



(17) 

where is the viscous stress in the (/(-direction as explained 
above. 

The last conservation equation to consider is that for energy 
per unit area within the SL. This is composed of many pieces, 
the first being kinetic energy, 



E KE = ^(v 2 g + vl)y. 
The gravitational potential energy is 



^grav 



-L 



4 ges kTy 

PgeSZdz= . 

o 5 g eS - g md fim p 



The internal energy of the gas is 



■-gas ■ 



It 

2 J 



/o V m p 
and the radiative energy is 



pkT , 6 kTy 
■dz=-- 



5 fim p ' 



(18) 



(19) 



(20) 



^rad 



-/ 

Jo 



aT 4 dz = 



12 g md kTy 



5 geff-gmd V*n p 

The total radially integrated energy density is then E XoK - 
Eg m v + E g . ds + E m d. With this we express energy conservation 
as 



(21) 
Eke+ 



c) 



1 d 



— (2irRE tot cos9)+- — 



dt 



Rd9 



2ttRv(,cos9 E tot + / Pdz 



o 

■■-2ttRF cos 9, 



(22) 

where the extra term added to the advected energy density 
is the radially integrated work that is done as fluid elements 
move higher in latitude. This expression can be simplified to 
give 



dE tot d 

Rcos9 ^r + d9 



ve cos 9 



£ tot + J Pdz 



-RFcost 



(23) 



Equations (9), (1 1), (17), and (23) must now be solved to find 
the properties of the SL. 



3.2. Taking the Steady State Limit 

As in IS99 we solve the SL equations in the steady state 
limit (d/dt = 0) which greatly simplifies the solutions be- 
cause 9 is the only independent variable. To be a valid ap- 
proximation, we require a sink for the accreted material at 
some latitude above the equator. The material near the equa- 
tor moves in latitude rather quickly with v s ~ 10 5 - 10 6 cm s _1 
(as we show in §4) and thus passes through the SL on a 
timescale of tens of seconds. Above this region the solutions 
show a long, deep minimum in vg for most of the star where 
vg ~ 10 4 cm s" 1 . This means that fluid still moving higher in 
latitude can only reach the pole on a timescale of ~ 10 4 s. At 
this point the fluid will have lost its azimuthal velocity, dissi- 
pated its high internal energy, and have an entropy set by the 
underlying WD luminosity. This means that long before the 
fluid has time to reach the pole it will mix with the slowly 
downward advecting fluid of the WD surface and thus be part 
of the WD. Recent theoretical calculations (Townsley & Bild- 
sten 2003) and observations (Sion 1999) show that accreting 
WDs have effective temperatures T e ff w (1 -3) x 10 4 K. At 
a latitude at which the SL has reached a similar T e ff and has 
« we therefore set a sink where we stop our integrations. 

The steady state assumption implies that continuity is set 
by the accretion rate, so that integrating equation (9) results 
in 



2 



2nRyvg cos 9. 



(24) 



We use y vg cos 9 = constant and the steady state approximation 
to simplify the other conservation equations. Equation (11) 
becomes 



dvg 4 d 
yvgcos9—+-cos9— 
d9 5 d9 



geff kTy 



L£eff 

= —Rt$ cos 
whereas equation (17) is rewritten 



"grad pm p 

'-yv£sin0. 



7?T0 cos 6 +yv$v r f, sin 9. 
Equation (23) is simplified to 



yvg cos9—f 
d9 



(25) 



(26) 



yvg cos9- 



vl + vl 



2 7g eff +3g rad kT 

5 geff-grad M m P 

= -RFcos9. 



(27) 

Rewriting equations (25), (26), and (27) in terms of the three 
dependent variables, vg, v^, and T we get, 



,dvg 



4 
— i 
5 



vi — - + -vi cos9— 



gets kT/(jjm p ) 



geff-grad VgCOS0 
---Fg-F cen , 



VgVtf, 

7 dvg 2 d 
V °^9 + 5 Ve d9 



■ — —Fa, + F„ 



dv q 
d9 

7g e ff+3g rad kT 



(28a) 
(28b) 



geff-grad _ 

=F^-F cen -E ai , (28c) 

where the terms appearing on the right sides are: the frictional 
dissipations in the <p- and ^-directions, 



Fg 



R 2 
a-v t 



vl + vl, 



R 2 R 



(29) 



(30) 
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where h = P/(pg e s) = kT /[/j,m p (g &s - g mA )] is the pressure 
scale height, the centrifugal term, 



F cen = VgV 



2 tan 6>, 



and the radiative dissipation 



RF M Edd 

£rad = = RgmdVe —T~ COS 6, 

y M 



(31) 



(32) 



where M^id = 4-tvRc/k is the Eddington limit (with no rota- 
tion). In Appendix B we show how to solve equations (28 a)- 
(28c) for each of the derivatives of the dependent variables to 
assist numerical integration. 

4. SOLVING FOR THE SPREADING LAYER PROPERTIES 

We now solve equations (28a)-(28c) by integrating from 
the equator toward the pole and halting integration once the 
sink has been reached. There are three boundary conditions at 
the equator for each of the dependent variables: vg o, v^o, and 
To. Furthermore, we need an initial angle at which to begin 
our integrations, 9 . In §4.1 we consider the range of values 
appropriate for each of these variables, and show that because 
the accretion rates are so much less than the Eddington limit, 
the integrations naturally attract toward solutions where local 
viscous dissipation is approximately balanced by local radia- 
tive cooling. In §4.2 we analytically consider the differential 
equations to gain some intuition before considering the nu- 
merical results, including showing how #sl changes with M 
and a. In §4.3 we present the full numerical solutions with 
descriptions of the important SL properties. 

4. 1 . Boundary Conditions at the Equator 

We begin by assuming that the accretion disk can be ap- 
proximated by a standard, thin Shakura-Sunyaev disk. Using 
the disk height from Shakura & Sunyaev (1973; in the limit 
of gas pressure much greater than radiation pressure and using 
Kramer's opacity) and dividing by R we find the disk subtends 
an angle at the WD surface of 

^ sk =1.5xlO- 2 ^ 10 MrM[ 3/8 < 8 , (33) 
where ayisk is the viscosity parameter for the accretion 
disk expected when the disk is actively accreting, Mig = 
M/(10 18 g s" 1 ), Mi = M/Mq, R 9 = R/(10 9 cm), and we set 
the factor [\-{r/R)~ l l 2 ] « 1. We have chosen a rather high ac- 
cretion rate for normalization because this is the regime where 
the SL will have the most observational impact. The material 
which comes in from the disk and subsequently spreads over 
the star must then start at an angle 6q < #disk- We therefore 
use a starting angle for integrations of 9o « 10~ 3 . In cases 
where the spreading angle is comparable or less than #disk we 
say that spreading is negligible and most of the dissipated en- 
ergy is radiated back into the disk. For the initial azimuthal 
velocity we use v^o = 0.99vk so that the accreted material is 
moving in nearly Keplerian orbits when it reaches the star. 
This leaves two parameters to determine, v^o an d ?b- 

The initial temperature of the SL is set using the disk's 
midplane temperature close to the WD. This is found from 
Shakura & Sunyaev (1973) to be 



(34) 



Finally, we must set the value of v#,o- To see that there is a pre- 
ferred value for this, consider setting v^o, To, and O using the 
description above. We then perform a number of integrations, 




0.001 



0.01 



FIG. 2. — A number of sample integrations used to constrain the boundary 
condition at the equator. For all of the integrations we fix the values of 9q = 
1CT 3 , v^o = 0.99v K , a = lfr 3 , a disk = 0.1, andM = 3 x 10 17 g s" 1 , for a white 
dwarf with M = 0.6Mq and fi = 9x 10 8 cm. Tq is fixed using equations (34). 
Only vgo is slightly different each time. Looking at the top panel it is clear 
that there is one preferred solution for vg (dotted line). The bottom panel 
shows the ratio E rd<i /F^ for the same set of integrations. Comparing the two 
panels, we see that this preferred solution corresponds to E ra d m F^. 



keeping all boundary values fixed while slightly varying v^o 
each time. The results of these integrations are shown in Fig- 
ure 2. In the top panel we see that the vg is naturally attracted 
toward a specific solution. 

To find out the physical meaning behind this preferred so- 
lution, we consider these same integrations, instead following 
the ratio E^/F^, in the bottom panel of Figure 2. This ratio 
reflects how efficiently the system is locally radiating the en- 
ergy dissipated from friction. If E md < F<p then there must be 
advection because there is more local viscous heating than can 
be radiated away. On the other hand, if E a & > F^ the mate- 
rial loses internal energy because it is radiating so efficiently. 
(We can ignore Fg because it is so much less than F^.) Com- 
paring the top and bottom panels, we see that the asymptotic 
solution found for vg appears to correspond to the case when 
initially E rid = F^, which is when all of the viscous dissipa- 
tion is being radiated away locally. To investigate further, we 
take a closer look at E^/F^ in the top panel of Figure 3. This 
shows that in fact the attracting solution is E md < which 
means that a small amount of advection must take place. In 
the bottom panel of Figure 3 we see that the actual ratio de- 
pends on the accretion rate and that the higher the accretion 
rate is, the more advection takes place. With this information 
in hand, we then use £ ra{ j < F$ to get a closed set of equations 
so that the SL can be solved uniquely for a given M and a (see 
Appendix C for additional details). 

This boundary condition is different than the NS case stud- 
ied by IS99 where it was found that the solutions are attracted 
toward the Eddington limit (i.e., g e ff = g ra( j) and there is a great 
deal of advection. The differences between these results high- 
lights the two limits of SL solutions. In the case of NSs, the 
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FIG. 3. — The top panel shows the same integrations as in Figure 2, but 
focusing near the equator. This shows that the actual asymptotic solution has 
EraA < F^,, so that there is naturally a small amount of advection taking place. 
The bottom panel shows that solutions for accretion rates of 10' 7 , 3 X 10 17 , 

5 X 10 17 , and 10 18 g s _1 (from top to bottom) with all other values the same as 
in Figure (2). All solutions show £ ra d < F^ and furthermore, this inequality 
increases with higher accretion rate. As the accretion rate increases toward 
the Eddington limit we expect a great deal of advection as found by Inogamov 

6 Sunyaev (1999) for the neutron star case. 



local accretion rate through the disk is much closer to the Ed- 
dington rate making advection important. The low M limit 
is the case of WDs considered here where radiation pressure 
is negligible so that little advection takes place. A similar 
regime likely exists for NSs at low accretion rates as well. 
In Figure 4 we show the initial latitudinal velocity, v^o, us- 
ing the condition E md = along with g t <s = g m &- At low 
accretion rates the two conditions are fairly different and the 
WD SL solutions are shown to have little advection and differ 
greatly from the NS case. As the accretion rate increases, the 
two conditions asymptote together and the WD solutions be- 
come more like the NS case. For boundary values at the right 
of the maximum in vg,o much more advection is seen in the 
solutions. 

Given the prescription described above there is still an addi- 
tional constraint for the values M and a to be consistent with 
our plane parallel model which assumes that h <C R- Substitut- 
ing the boundary condition that E md = F^ into h < R provides 
the constraint 



< 10 



2-3/5,^/10^/4^/4 



(35) 



This shows that at high M, or when there is a lot of viscous 
dissipation because of a high a, it is possible that the SL will 
expand so that h w R. This case is interesting for some high 
M systems (such as symbiotic binaries which can have M w 
10 20 g s" 1 ), so although outside the scope of this paper, such 
regions of parameter space should be studied in the future. 

4.2. Analytic Estimates 
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FIG. 4. — The initial latitudinal velocity, vg Q , versus accretion rate, using 
the boundary condition that £ ra( j = (solid lines). The dashed line denotes 
the Eddington limit g e f t = g a i and the dotted lines pass through the maximum 
value vgfi can attain for a given viscosity, a. Solutions taken from the left 
of the dotted line show little advection and a strong attraction toward the 
E m d = F^. Solutions using v e from the right of the dotted line show much 
more advection, similar to the case of neutron stars investigated by Inogamov 
& Sunyaev (1999). 



As a first attempt at approximating the spreading angle, 6*sl, 
one might guess that it is given by geostrophic balance (when 
the pressure gradient in the latitudinal direction is balanced by 
the Coriolis force pushing back toward the equator). This can 
be found using conservation of momentum in the (^-direction, 
equation (25). Approximating v <C v^, g rd d < g e ff. and 9 <C 
1, this results in 



d kTy 
d6 [im, p 



This then gives #sl approximately as 



0: 



SL 



kT 

fim p 



1/2 



(36) 



(37) 



This may provide the correct estimate when friction is neg- 
ligible, but it gives the incorrect scalings for the case we are 
considering here where friction is the primary agent for slow- 
ing the flow. 

Using this insight, we now derive approximate solutions, 
instead assuming that conservation of 0-momentum is con- 
trolling the flow. This will illuminate the general features of 
the SL and provide power law scalings consistent with numer- 
ical integrations. We first must estimate the size and scaling 
of v$. Setting E m & = F^ near the equator we find 

JW fl -^»a-v^/vg + v£»a-V (38) 

Solving equation (38) for vg gives the approximate latitudinal 
velocity in the SL. Using equation (34) to set the temperature 
we find 

v e , SL = 2.1 x 10 6 cm s-'a^af M} g /4 Aff /8 R~ V \ (39) 
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where a 3 = a/10~ 3 . 

To understand the effect of friction on setting 6*sl we use 
conservation of momentum in the 0-direction, equation (28b), 
which near the equator is 

dvd, _ R , 



Vgv^- 

We can then estimate 



dO 



h V0V4, 
R av\ ' 



(40) 



(41) 



This gives the simple explanation that 6*sl is set by two fac- 
tors: (1) the pressure gradient in the latitudinal direction (the 
first term), and (2) the competition between (^-momentum 
fluxing the in the ^-direction, vgv^, and dissipation of <jh 
momentum due to viscous stress, avj. This scaling is dif- 
ferent than the results of IS99 who found that the spreading 
flow adjusts itself so as to radiate everywhere at close to the 
local Eddington rate, and thus their #sl depends only on the 
ratio M/Mem- At much less than the Eddington rate, equation 
(41) shows that #sl is proportional to h and therefore highly 
dependent on the initial properties of the material when it first 
reaches the WD equator. We can find the functional relation- 
ship between and 9 by assuming that vg and T are constant 
and integrating equation (40). Setting g eff = GM/R 2 -v^/R 
and using 8 w this results in 



6 = 



1 kT Vg 

a nm p v\ v K 



x ( i* - J5L + log vg5gj _ l 0g V^Wj ) (42 ) 

\ v v <t>,0 V 1 + V 0/ V K V 1 + V W V K/ 

Figure 5 shows this result plotted alongside the numerical so- 
lutions found in the next section. The analytic estimate shows 
only qualitative agreement, which is not surprising given the 
crudeness of the estimates made, namely that both vg and T 
are constant. We estimate the spreading angle, #sl, by assum- 
ing that w 0.1 Vk at the end of the spreading. In the limit of 
v<t>,o ~ v k we find 



1 kT Vg ( 1 18/11 
&sl = r — 9 + -log- — 

a nm p v^ v K V 2 1 ~ v <P,o/vk 



(43) 



This result has a singularity at v^o = vk (which it must be- 
cause material moving at the Keplerian velocity will not settle 
onto the WD), but even for values v^o < vk this last term is 
only of order a few and has little effect on this angle estimate. 
Substituting v g and T by using equations (39) and (34), and 
setting v^ o = 0.99v K we find 

n <-> / irv-2 3/10 -1/2 .>1 1/20, .-5/8 „-l/8 

6> SL = 2.4 x 10 z aJ isk a 3 1 M l& < M l ' R g 1 . (44) 

The scalings predicted by equation (44) are consistent with 
those found performing the numerical integrations (as we de- 
scribe in the following section). This derivation shows that 
the change in v^ and the value of #sl are mainly dictated by 
viscous dissipation. Using (d/dt)(yv<f,) = -apv^, we estimate 
a timescale for spreading of 

f SL = h/iavj,) w 50 s a^aj 1 ^ 10 ^ 4 ^ 4 . (45) 

This timescale is extremely short and mainly due to the small 
amount of material which is in the SL at any given time. This 
implies that if accretion decreases dramatically (such as fol- 
lowing a dwarf nova outburst) the hot belt around the equator 
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FIG. 5. — A comparison between the numerical integrations (solid lines) 
with the analytic estimate given by equation (42) (dotted lines). Accretion 
rates of M = 10 17 , 10 18 , and 10 19 g s~' are shown (from left to right), all 
with a = 1(T 3 , o disk = 0.1, M = O.6M , and R = 9 X 1() 8 cm. Though quan- 
titatively different, the two results are qualitatively similar and both predict 
similar spreading angles at = O.Ivk ~ 3 X 10 7 cm s _1 . At higher accretion 
rates the solutions differ more because the analytic estimate ignores radiation 
pressure effects. 



will spread fairly quickly over the WD surface and not be vis- 
ible days or even hours later. 

Since at low accretion rates #sl ^ #disk, the spreading ma- 
terial will most likely be obscured by the disk, making direct 
observations of the SL difficult. Setting 9sl = Odisk we solve 
for the critical M at which the SL will extend above the disk, 



M crit > 3.7 x 10 17 gs 



- 1 a 3/ V / V 5/ V /8 



(46) 



The large value of this accretion rate implies that the effects 
of the SL may only be seen in a select number of systems 
where a high M occurs such as dwarf novae, symbiotic bi- 
naries, and supersoft sources. Lower accretion rate sources 
must still have dissipation occuring at the interface between 
the disk and WD surface, but the radiation created in these 
cases is most likely absorbed and re-radiated by the accretion 
disk. We discuss all of these issues in more detail in §5. 

Using the results from above we find that the local flux in 
the SL is approximately, 



Fsl = 1 x 10 erg cm" 



-2 -1 -3/10 1/2.V9/20, , 13/8^-23/8, 
S "disk "3 M A M l R 9 (47) 



Multiplying by the area of the spreading zone which is 
4ttR 2 sin 6>sl ~ 47t/? 2 6>sl provides the luminosity. Comparing 
equations (44) and (47) we find that L oc M and independent 
of both a and oyisk as it must be since these are free param- 
eters. The value of the luminosity is consistent with standard 
energy conservation arguments, L = GMM/(2R). This means 
that luminosity alone will not tell us whether a SL is present. 
A better test would be a study of the emitting area and T e ff. 
The emitting area scales like 0$l which is expressed above. 
The T e ff is given by the flux so that from equation (47) we 
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FIG. 6. — The velocities and temperature of the spreading layer for a white 
dwarf with M = O.6M and R = 9 X 10 8 cm, with a = 1(T 3 and a disk = 0.1. 
From left to right the accretion rates are 10 17 , 10 18 , 10 19 and 3 X 10 19 g s -1 . 



find 

7k, SL = 2 x 10 s K a-fi 4 \fMTM\ 3/3 W y32 - (48) 

This value is not so different than what would have been ap- 
proximated from BL arguments (for example see Frank, King 
& Raine 2002). This is not surprising since even if the SL in- 
creases the emitting area by a factor of 10, the effective tem- 
perature would only change by a factor of 1 .8. The best test to 
see if spreading is taking place is to follow how T e ff changes 
with M. The SL model shows that T e g will change much more 
gradually than if the radiating area on the WD is fixed (which 
would predict T e ff oc M 1 / 4 ). 

4.3. Numerical Solutions 

We now present the full numerical solutions of equations 
(28a)-(28c). These confirm the general estimates made in 
the previous section and also provide additional details that 
cannot be investigated analytically. For these integrations we 
use the fiducial values of M = 0.6M Q , R = 9 x 10 8 cm, and 
adisk = 0.1. The composition of the accreting material is as- 
sumed to be solar so that [i = 0.62 and k = 0.34 cm 2 g _1 . 

First we set our boundary value for To using the temperature 
of equation (34), and then we set v^o by using the condition 
that E m d < F$. The integrations are then performed on the 
three coupled differential equations as derived in Appendix 
B from equations (28a)-(28c). In Figures 6 and 7 we plot 
the properties of the SL for a = 10" 3 and a range of accre- 
tion rates. At low accretion rates (M w 10 17 - 10 18 gs" 1 ) we 
find the numerical integrations result in values close to what 
was approximated analytically in the previous section. At the 
higher accretion rates (M > 10 19 gs" 1 ) more advection takes 
place, resulting in solutions which look similar to what IS99 
found for neutron stars. Most notably, the temperature pro- 
files show large peaks at high latitudes due to energy being 
advected away from the equator and deposited higher up on 
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FIG. 7. — The effective temperature, optical depth, and pressure scale height 
of the spreading layer, using the same parameters as in Figure 6. The dotted 
line shows T c ff = 3 X 10 4 K, a fiducial temperature for the underlying accret- 
ing white dwarf. 
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FIG. 8. — The free-free opacity, Kff, calculated using the properties of the 
spreading layer solutions, using the same parameters as in Figure 6. This can 
be compared to the electron scattering opacity, k C s = 0.34 cm 2 g _1 (dotted 
line), to see that the opacity is correctly approximated for the majority of the 
spreading layer, especially in the regions which provide the majority of the 
flux. 



the star. In Figure 7 we denote the fiducial T e s expected for the 
underlying accreting WD. Comparing this and Figure 6 shows 
that once the SL has come into thermal balance with the un- 
derlying WD, the SL is also rotating with the WD (v<f, « 0). 
The optical depth is shown to always be much greater than 
unity, so that our assumptions of a radiation pressure acting in 
our one-zone envelope and a true T e g are valid. At high accre- 
tion rates the scale height becomes comparable to R, show- 
ing the limit of the plane parallel model. Our viscous stress 
parametrization becomes questionable at low accretion rates 
(M < 10 17 g s" 1 ) when the latitudinal shear is much greater 
than the shear between the SL and underlying WD (because 
R9sl < h). Qualitatively, this does not change the important 
result that most or all of the SL dissipation happens within the 
width of the accretion disk at these accretion rates. Figure 8 
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FIG. 9. — The velocities and temperature of the spreading layer for a white 
dwarf with M = 0.6M Q and R = 9 X 10 8 cm, with M = 10 18 g s _1 and «disk = 
0.1. From left to right the values of a are 3 X lfT 3 , 1(T 3 , 3 X lfT 4 , and Kr 4 . 



FIG. 10. — The effective temperature, optical depth, and pressure scale 
height of the spreading layer, using the same parameters as in Figure 9. 
Viscous stress decreases from left to right. The dotted line shows T e ff = 
3 X 10 4 K, a fiducial temperature for the underlying accreting white dwarf. 



shows the free-free opacity, calculated from the solutions that 
are found using a constant electron scattering opacity. For all 
solutions, the free-free opacity is small in the hot portions of 
the SL. This gives some confidence that T e ff over this region 
is accurate, because the majority of the flux is radiated here. 
Unfortunately, our opacity assumption always breaks down in 
regions where the SL is cool which makes any spectrum cal- 
culation based on this model somewhat dubious. In Figures 
9 and 10 we show that the properties of the SL can vary de- 
pending on the value of the viscosity. 

Figure 1 1 summarizes the key feature of the SL which is 
how the spreading angle (and thus the radiating area) changes 
with accretion rate. At low accretion rates the scaling is con- 
sistent with what was found analytically, equation (44). Fig- 
ure 1 1 reinforces the result that high accretion rates are needed 
for any hope of appreciable spreading to take place. It is im- 
portant to note that even when the accretion disk covers the 
main portion of the SL, the SL may still have some observa- 
tional impact depending on the viewing angle of the WD and 
how far the contrast in r e ff extends (see Figures 7 and 10). 

5. INITIAL COMPARISONS TO OBSERVATIONS 

From the analysis in §4, it is qualitatively clear that the SL 
is best observed in high accretion rate systems so that it is not 
obscured by the accretion disk or its winds. This suggests that 
the best candidates to show a SL are dwarf novae, symbiotic 
binaries, and supersoft sources. We now review the current 
observations of these accreting WDs in the EUV and soft X- 
rays. 

5.1. Dwarf Nova Outbursts 

There have been a number of observations of dwarf no- 
vae in outburst or superoutburst using the Extreme Ultravi- 
olet Explorer {EWE), including those of SS Cyg (Mauche, 
Raymond & Mattei 1995; Wheatley, Mauche & Mattei 2003, 
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FIG. 11. — The spreading angle, #sl. versus accretion rate for a white 
dwarf with M = 0.6Mq and R = 9 X 10 9 cm and three different viscous stress 
parameters. The spreading angle is defined as the angle at which = O.Ivk- 
At low accretion rates the scaling with M is fairly consistent with what is 
expected from an analytic approximation which gives 6*sl <x M 11 / 20 - This 
also shows that high accretion rates are needed for the spreading layer to 
extend beyond the angle made by a Shakura-Sunyaev disk, equation (33) 
(dashed line). 



10 



hereafter WMM03), U Gem (Long et al. 1996), VW Hyi 
(Mauche 1996) and OY Car (Mauche & Raymond 2000). 
Mauche, Mattei & Bateson (2001) provides a review of many 
of these objects. 

WMM03 provide an especially interesting observation of 
an outburst of SS Cygni studied in the optical by The Amer- 
ican Association of Variable Star Observers (AAVSO), in the 
X-rays by the Rossi X-ray Timing Explorer (RXTE), and in the 
EUV using EUVE. The observations show a complete optical 
light curve, following the burst over ~ 12 days. The X-rays 
show two small bursts of activity, the first coming ~ 1 .5 days 
after the initial rise in the optical and the second coming near 
the end of the outburst. In between there is emission in the 
EUV. The correlations between these three bands are consis- 
tent with what is expected by the standard dwarf nova picture. 
The outburst begins once the disk is in the high state, start- 
ing the rise in optical, and then proceeds from the outside of 
the disk in. Once material begins accreting onto the WD it is 
initially optically thin, and thus radiates in the X-rays at a tem- 
perature given by the virial temperature of the accreting gas 
(Pringle & Savonije 1979; Tylenda 1981; Narayan & Popham 
1993). This provides the initial <~ 1.5 days of X-rays. Once 
the material becomes optically thick, it radiates in the EUV 
for most of the remainder of the burst (Pringle 1977; Tylenda 
1977; Popham & Narayan 1995). Finally, at the end, the ma- 
terial becomes optically thin one last time, radiating a little 
more in the X-rays before the outburst ceases. 

Qualitatively, a similar picture is also expected for a SL, 
even though the calculations we perform above only apply 
when the SL is optically thick. One might expect the SL 
to be important because of the large inferred accretion rates 
at the peak of the outburst (~ 10 18 g s -1 ). As discussed in 
§4.2, evidence that spreading is occuring requires studying 
the change in effective temperature with accretion rate. If the 
radiating area does not change with accretion rate then one 
expects T e ff oc M 1 / 4 . On the other hand, since the area in- 
creases with M for the SL we find T e ff oc M 9 / 80 , equation (48). 
Figure 12 shows the count rate (which is proportional to lu- 
minosity) versus the hardness ratio (related to T e ff) taken from 
this observation by WMM03. Fit are the two power laws ex- 
pected from a SL and a constant radiating area. Even though 
this does not show conclusively that the SL exists, this shows 
in principle the kind of test which needs to be done as data 
improves. Other observations of SS Cyg (such as Mauche, 
Raymond & Mattei 1995), also report a fairly constant hard- 
ness ratio during outburst even though the count rate changes 
by two orders of magnitude. 

Another suggestive result by Mauche, Raymond & Mattei 
(1995) is the count rate in the EUV as compared with typical 
count rates in the optical for similar outbursts. If the former 
is assumed to represent the region up close to the WD and 
the latter the accretion disk, these observations suggest that 
the ratio of the two luminosities is Lbl/Lh^ < 0.07. This 
is in stark contrast from LBL/^disk ~ 1 as expected from BL 
models. It is possible that because the SL only extends a small 
distance above the disk, much of the radiation from the SL 
is obscured by, or contained within, the accretion disk or its 
winds. 

A much more complicated picture is provided by the su- 
peroutburst of the edge-on binary OY Car (Mauche & Ray- 
mond 2000). By comparing optical and EUV lightcurves dur- 
ing outburst they find that while the optical is eclipsed, the 
EUV shows no sign of being eclipsed. This calls into question 
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FIG. 12. — Hardness ratio versus count rate as measured during a dwarf 
nova outburst of SS Cyg by Wheatley, Mauche & Mattei (2003). The hard- 
ness ratio is defined as the ratio of counts between the the 72 — 95 A band 
(excluding 76 - 80 A) to the 95 - 130 A band. Plotted over the data are the 
relations expected from the spreading layer (dotted line) and a model which 
has constant radiating area (dashed line). In each case the area is assumed 
to radiate like a blackbody and the Wein limit is used. The scaling between 
count rate and accretion rate is treated as a free parameter used to maximize 
the fit in each case. The spreading layer gives a shallower power law, as ex- 
pected because its radiating area increases with accretion rate, but it provides 
only a slightly better fit to the data. 



the assumption that the EUV is associated with a region near 
the WD and suggests that perhaps the EUV is from a much 
more extended region than previously thought. Good fits are 
made to the outburst spectrum by assuming that the radiation 
from both the BL and accretion disk are scattered into the ob- 
server's line of sight by an ionized accretion disk wind. These 
additional complications make it difficult to simply compare 
our SL with this or similar observation. 

In a study of U Gem using the Far Ultraviolet Spectroscopic 
Explorer (FUSE), Froning et al. (2001) model the spectra 
both during outburst and in the decline afterwards in the far- 
ultraviolet. They find that the majority of the spectra during 
outburst is consistent with a steady state accretion disk with 
M = 4.6 x 10 17 g s" 1 . This is close to the critical accretion rate 
found in equation (46) so that we might guess a SL is present. 
The accretion disk model underpredicts the observed spectra 
at short wavelengths (< 950A), which is a sign of a missing 
hot component, perhaps the SL. Though suggestive, there is 
still more work to be done to determine whether a SL has 
been seen. It is not clear whether accretion disk model im- 
provements would explain this difference. In the same study 
by Froning et al. (2001), they also look at U Gem a few days 
after the outburst peak when the luminosity is declining. They 
find that the spectra are very well fit by a WD model. They 
also attempt a fit which includes an additional hot belt on the 
WD, but they find that the improvement of the fit is negligi- 
ble. This is long after the timescale we find for spreading, 
equation (45), so we do not expect to see a SL in this case. 
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In an observation of VW Hyi, Sion et al. (1996) use the 
Hubble Faint Object Spectrograph (FOS) to look at spectra 
a day after a dwarf nova outburst and ten days after a super 
outburst. In each case they compare spectral fits of a single 
WD model with the fits of a WD plus a hot belt, and for both 
observations they find that the latter model provides a better 
fit. From the solutions found for the SL in §4 it is surprising 
that the hot belt is still present on the WD this many days after 
outburst when the SL should have already spread over the star 
due to friction. Other observations by the same group (Huang 
et al. 1996a, 1996b) show similar, but less conclusive results. 

5.2. Symbiotic Binaries and Supersoft Sources 

Another promising candidate for showing a SL are symbi- 
otic binaries. These are WDs which are accreting from the 
winds of red giants with M = 10 17 - 10 20 g s" 1 , and they typ- 
ically show a hot component which may or may not be com- 
ing from near the WD. Miirset, Nussbaumer, Schmid & Vogel 
(1991) review the characteristics of this hot component for 
a large number of systems, finding luminosities in the range 
0.3L© < L < 3.7 x 10 4 Lq and effective temperatures around 
5.5 x 10 4 K < r eff < 2.1 x 10 5 K. From these two quanti- 
ties the radius associated with the emitting region can be es- 
timated. This is typically ~ 10 10 cm, much larger than ex- 
pected if the radiation where coming directly from the WD 
surface, which suggests the presence of a hot ionized nebula 
surrounding the WD. There are three systems listed in Miirset 
et al. (1991) which have features qualitatively different than 
most other symbiotic binaries; these are EG And, AG Dra, 
and CH Cyg. All three of this systems are underluminous in 
comparison with the other symbiotic binaries which implies 
a much smaller radius, closer to what is expected for WD. 
A May 1988 observation of CH Cyg resulted in an inferred 
radius of 1.4 x 10 8 cm which is even too small to be a WD. 
Though Miirset et al. (1991) argue this could be due to an 
incorrect distance measurement, it could also be a sign of the 
SL. Two of these systems show flickering in the optical (EG 
And and CH Cyg; Sokoloski, Bildsten & Ho 2001), which is a 
trait uncommon to symbiotics, but common to accreting CVs. 
This may point to disk accretion being present and thus also a 
SL. 

Finally, supersoft sources are WD binaries with M <~ 10 18 - 
10 20 g s -1 . At these high accretion rates accreted hydrogen 
on the WD surface is believed to be burning steadily (van 
den Heuvel et al. 1992; Heise, van Teeseling & Kahabka 
1994) and this provides the large luminosity seem from these 
sources (~ 10 37 - 10 38 erg s" 1 ). The SL would appear at a sig- 
nificantly different T e ff from the WD surface because of the 
difference in emitting area. The problem with these systems 
is that the burning causes such high luminosities that the SL 
would be dwarfed in comparison. Current observations of su- 
persoft sources do not have enough signal to show a SL and 
would require much more detailed models to correctly deter- 
mine whether a SL is present. 

6. DISCUSSION AND CONCLUSIONS 

The SL model, first conceived by IS99, provides a new 
way of understanding the properties of newly accreted ma- 
terial and its interaction with the stellar surface. In this model 
the latitudinal direction is used as the independent variable, 
as opposed to standard BL models which follow how the disk 
changes near the star along the radial direction. This allows 
an investigation of the radiating area of the hot belt which 
forms near the equator, and it describes the properties of the 



accreting material when it first comes into hydrostatic bal- 
ance on the star. We have investigated the solutions of the 
SL model when applied to WDs, accreting in the range of 
M ~ 10 17 - 10 19 g s" 1 . We find that the integrations for the 
spreading flow naturally attract toward solutions in which ini- 
tially local radiative cooling balances local viscous dissipa- 
tion. As the accretion rate gets higher, we find increasingly 
more advection in the flow, similar to the solutions for NSs 
investigated by IS99. The spreading angle on WDs is found to 
typically be 6>sl ~ 0.01 -0. 1, depending on both the accretion 
rate and the viscosity. The SL has an effective temperatures 
in the range of (2- 5) x 10 5 K with a pressure scale height of 
10 7 -10 8 cm. 

To be clearly observed in actively accreting systems, the 
SL must extend to a large enough angle to be seen above the 
accretion disk. This requires the accretion rate to be high, 
M > 10 18 g s" 1 , assuming that the accretion disk is similar to 
a thin Shakura-Sunyaev disk. The best candidates to show the 
SL are dwarf novae in outburst, supersoft sources, and sym- 
biotic binaries. Current observations of the dwarf novae may 
be the best opportunity to see if spreading is actually occur- 
ring. The scaling of effective temperature with accretion rate 
is much weaker in the case of a SL because of the change in 
radiating area. At lower accretion rates, M < 10 18 g s" 1 , the 
SL will not spread far from the equator and most of the dis- 
sipated energy will radiate back into the accretion disk. Al- 
though this may make comparison between theory and obser- 
vations difficult, because the SL cannot be seen directly, it 
may still have important consequences for the accretion disk 
structure and spectra. 

The short timescale for spreading, ?sl ~ 50 s, suggests that 
a hot belt around the WD equator will be difficult to see after 
an outburst because the hot material should spread over the 
WD surface quickly. This also means that mixing between 
freshly accreted material and much deeper layers appears dif- 
ficult on accreting WDs. Rosner et al. (2001) and Alexakis 
et al. (2003) have argued that the mixing of recently accreted 
H/He with the underlying layer of C/O due to rotational shear- 
ing explains observations of CNO nuclei in the ejecta of nova. 
To have ?sl of order the accretion time (~ 10 yr) requires 
a < 10~ n , so that the viscosity must be even less than what is 
estimated from microphysics (see Appendix A). This makes 
it difficult to imagine mixing via this mechanism all the way 
down to the C/O layers. 

The initial calculations by IS99 and in this paper describe 
some of the general features of the SL model, a new area of in- 
vestigation which will lead to further studies of how accreted 
material settles onto stars. This may include studying how 
differential rotation and angular momentum transfer affect the 
underlying stellar surface or the possibility of nonradial oscil- 
lations present in the SL. Spectral modeling of the SL, along 
comparisons with observations, may also help in identifying 
if and when spreading is present. 
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APPENDIX 

A. VISCOUS STRESS PARAMETRIZATION 
In the derivation of the spreading differential equations we choose to parametrize r in terms of a unitless constant, a, so that 



= apv 



(Al) 



We follow IS99 and estimate the values of a implied by possible ionic or radiative viscosities. In general the viscous stress is 
expected to be set by turbulence, so these estimates only provide a lower limit to a. The ion viscosity, z/,-, is related to the viscous 
stress by 



r = Vip 



dv 
dz 



v iP-, 



(A2) 



Using the ion viscosity of Spitzer (1965) we find Vi ss 2 x 10 6 cm 2 s" 1 (using typical values of the spreading layer solutions from 
§4). Equations (Al) and (A2) can be combined to find a ~ 10~ 10 , which shows that this contribution can be ignored. We next 
consider the radiative viscosity, u r , which is given by 



r = v r p r 



dv 



v r p 



(A3) 



where p r = aT 4 /c 2 is the radiation density. The radiative viscosity can be be written as v r = Ac, where A is the photon mean free 
path given by A = l/(pn). We then solve for a and find it can be written in terms of three ratios as 



a = 



Pr AC 

p hv 



1(T 



(A4) 



which is small due to p r j p <C 1 . Though larger than the former value for a, this quick check shows that neither of these 
mechanisms are large enough to provide the viscosity needed to slow down the spreading material. We therefore assume that 
there is a turbulent viscosity which results in higher values of a and leave this as a free parameter. In IS99, this form of the 
viscous stress, equation (Al), is contrasted with the form typically used for accretion disks (r = ac s h, where c s is the speed of 
sound; Shakura & Sunyaev 1973). The latter prescription for the viscous stress describes the friction between adjacent annuli of 
the layer and IS99 show that it is negligible as long as h/R <C 1 (so that most of the free energy is at the SL/WD boundary). 

B. REWRITING THE DIFFERENTIAL EQUATIONS FOR NUMERICAL INTEGRATION 

To solve equations (28a)-(28c) numerically we must solve for each of the derivatives in terms of the dependent variables. This 
involves expanding the messy derivatives in equations (28a) and (28c). From equation (28a) we need 



d 



ges kT/(pm p ) 



geff "grad VgCOsO 



geff 



1 d kT g eS kT/(pm p )dv g 



1 



kT/(pm p ) dg e « 



g s ff- g rd d V0 COS0 dd pm p geff-grad vjcOsO 

geff kT/(pm p ) d 



dO geff "grad Vg COS 6 d6 



(geff-grad) 2 Vg COS 6 d6 

The individual derivatives in this expression are given by 

dg e ff 2ve dvp 2v0 dv<j, 

~df~ ~Y~d7 ~R~dJ' 

and 

dgrad 



Cgeff-gradH 



g eff kT lipnip) 

geff-grad VQ COS 9 



tan 6. 



d9 



4 dT 1 dv e , 
T dO vg do 



(Bl) 



(B2) 



(B3) 



We can then rewrite equation (28a) as 



4 geff " 2geffgrad ~ 2grzdVg /R kT 



(geff-grad) 2 



prrip 



grad 



dv B + 8 

dO 5 (geff "grad) 2 



v v0 kT dv<f 



R pm p d9 

= -Fg-F ce , 



4 geff+ 3 geffgrad , 
5 (geff" grad) 2 
4geff _ 2g e ffgrad 



■Vg 



v - 



5 (geff" grad) 2 pm p 



d kT 
dO pnip 

■t<mO. 



(B4) 



From equation (28c), we must expand 



7g e ff+3g rad kT 



geff" grad p m p 



10geff-3(g e ff-grad) kT 



geff "grad 



pm p 



10geff-3(g ef f-g ra d) d kT 

geff -grad dO pm p 

lOgeff 



(geff- grad) 2 pm p d9 



10 kT dg e n 

geff - grad M m p d 9 

kT d 

(geff "grad )■ 



(B5) 
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Using equations (B2) and (B3) we can then rewrite equation (28c) as 



2 , / \ gettg™d + 2g rdA v 2 e /R kT 

. 9 (geff-grad) 2 iOtlp 



dve 8g rad veV0 fcr dv,/, 



dO feeff-grad) 2 /? <^ 



2 7 £eff + 36g e ffg rad - 3g r 2 ad 
5 (geff-grad) 2 



d kT 

d6 /j,m p 



r „ p 4g ef fg rad fcr 

: i> - -Fcen - -Erad + : — ve tan 9 . 



(geff-grad)'" Vm p 



Equations (B4), (28b), and (B6) are of the form 



This can be inverted to find, 



dv dvj, d kT 
dv dv do fj,m p 

dvg dv,!, d kT 

B i h B 2 h D 3 = C3 . 

d6 d6 dO fintp 



dv _ B 3 Ci -A i C i + {A i B 2 -A 2 B i )C 2 /{ve V< p) 



AiB 3 -A 3 Bi 



dv,/, _ C 2 



d9 vgv^ 

d kT _A 1 C 3 -B l C l +(A 2 B l -A l B 2 )C 2 /(vgv 4> ) 



dO fjm p 

which can now be integrated numerically. 



A^-AiBi 



(B6) 

(B7) 

(B8) 
(B9) 

(BIO) 
(Bll) 
(B12) 



C. BOUNDARY CONDITION FOR LOW M 

Using the differential equations derived above, we can investigate the boundary condition vg.o in more detail. This reinforces 
that at low M the solutions naturally attract toward the boundary condition of local viscous heating being balanced by local 
radiative cooling. In the limit of 6 « 0, g e ff-g ra d ~ geff ~ GM/R 2 -v 2 i /R, and » ^JkT / (jim p ), vg, the terms from above 
become 



Ai 


2 4 kT 

= v e , 

5 /im ; , 


A 2 


=0, 


A3 


4 

= 5 V * 


Bi 


= v 2 e 


B 2 


= 


B3 


14 
= y v e 


c, 


= -F„ 


c 2 




c 3 


= ^0 — Erad • 



4 ytr 



5 [im p ' 



In this limit, the derivative of vg, equation (BIO), becomes 

14 17 4 



dvg _ ~—Fg 
~d6 



(F,p - E md ) 35F e +l 0(i> - £ rad ) 



%kT/{nm p )- 



4 2 

5 v e 



26kT I (jim p ) 



(CI) 

(C2) 

(C3) 

(C4) 
(C5) 

(C6) 

(C7) 
(C8) 
(C9) 

(CIO) 



Since Fg <C E m d, F<j,, the second term in the numerator will dominate unless E md = F$. If F^ > E md , then dvg/d9 > and the 
solutions will be driven toward higher vg. From continuity, equation (24), this means that y must decrease. The radiative cooling 
is £ ra d oc g ra dVe oc y~ 2 so that it will increase, pushing the solutions toward local energy balance. If E ml $ > F^, then dvgjdQ < 
and once again local energy balance is reached for similar reasons. The only time dvg jdO will vary slowly is when E m & = F<$, and 
we see in this case dvgjdQ « 35Fg/[26kT /(pm p )], which is positive. This is consistent with initial slope of vg in Figures (2) and 
(6). More detailed analysis, including the effects of g rad , is needed to correctly predict the values of the slopes seen. 
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